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Abstract: A new numerical approach to solving the classic problem of electromagnetic 
(EM) scattering off a perfect electric object is studied with the objective of 
substantially reducing computation times. The method considered here is the in the 
frequency domain Method of Moments (MoM) formulation involving the use of a dyadic 
Green's function (GF). Traditionally, this GF is formulated in free space, as afforded by 
the equivalence principle. However, since the resultant equivalent sources generate a 
null filed inside the scatterer volume, the door is open for the inclusion of arbitrary fillers 
therein. 


We suggest the usage of balanced absorbers as fillers and using their Green’s 
function instead of the free space one. To this end, the solution of the essentially volumetric 
problem of the absorber is required as a preprocessing stage. Balanced absorbers have both 
electric and magnetic Weston-like or Perfect Matched Layers (PML) loss mechanisms. Many 
interactions between pairs of basic functions are then virtually eliminated. As a result, the 
MoM matrix, representing the GF, is significantly thinned. 


The cost of calibrating this modified GF using the volumetric representation of 
the absorber is investigated. The effort incurred in the pre-processing stage can be alleviated 
by choosing absorber configurations that apply to many problems with high degrees of 
symmetry, thin absorbing shells rather than volumetric scatterers or homogeneous 
absorbers that lend themselves to surface formulations. 


It is shown that this form of thinning has little effect on the accuracy. Moreover, most 
of the thinned elements need not be computed at all. 
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Electric- 
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Absorber 


Green's Conventional surface integral equation formulations employ 


Method f equivalent sources with null field in the internal region of the 

ULT scatterer. Therefore, this region can be filled, at least 

ον conceptually, with a balanced absorber with both electric and 
magnetic conductances. The Green's function of the balanced 

absorbers is then employed instead of the conventional free space 

function. 

2734 By reducing interactions between pairs of opposing basis 

Sens functions, this method enables an MoM matrix with fewer than 

Conclusions 2596 non-zero elements with little effect on the result. Moreover, 


the annulled elements need not be computed. 
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zm The balanced condition is then 


Conditioning σ/ε — ο. (2) 


Example 


πο, Eq. (2) is sufficient for a 100% absorption of a normally incident 
plane wave at a planar interface between the homogeneous and 
absorbing medium. 


The idea of balanced EM absorbers is not entirely new 
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Vise The concept of the balanced absorber dates back to Salisbury in 


ne the 1940's! and to Weston in the 1960's?. It was later 
ite elaborated into the seminal PML by Bérenger. 


seine Balanced absorbers can work with very high or very low 
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ου ə The unknown source J(r’) is supported by the surface Co. 


Electric- 
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Heres e [n view of the null total field inside Co, we embed a 


UE balanced absorber in it. 


Function 


οί e The Green's function of the absorber, denoted Gar, r^), 


ULL and a corresponding modified incident field are to be used. 
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Electric- The rationale of using G54 is in the isolation it provides 
Magnetic ; ] 

ος between basis functions located across the absorber. As a result, 


(6εη’5 


"μέν ο) a typical row in the MoM matrix will decay faster away from the 


Method fi : ; 
W diagonal, relative to the conventional free space matrix: 




























Thinning T T T T 
Shay,  -——————————— ! b i : 
Kastner, = ] : 
Weile S Red: Free space matrix 
E —— Blue: Modified matrix — — 
Introduction 5 | | : | | : | 
o Bs 1i τετ τες Ετος οκ 5:ν''''''΄'.'.« νυν»: uie eee ti ee COME E 
BEMA (D | 
τ | 
Conditioning τ I ee ae ene μον ee es αρ μαμα παρ, continuam η 
© ! 
Example ΞΕ i 
ος ον ————————P— Án ((—————————————— 
: o | 
Conclusions zz : : 
5 | μάς | [i atari’ leseinst cs shes sd less 
0 50 100 150 : 200 250 300 350 


PECs and impedance surfaces have been used as fillers before^. 


^Boag, A. and V. Lomakin, 'Generalized equivalence integral equations," 
IEEE Antennas Wireless Propag. Lett, Vol. 11, 1568-1571, 2012. 
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Formulating an integral equation with the Balanced EM 
Absorber 


Reall: The balanced absorber is characterized by (2): 


c/e—oc"[y. 


The problem can then be formulated by the following IE 


] Gear. rar = ο ο (A. (3) 


Co C4 


Here, J™°(r’) is an equivalent source distribution that has 
generated the original incident field under free space conditions. 
It is now used to generate the modified incident field in the 
presence of the absorber, seen in the right hand side of (3). It 
may pay to choose C' as a straight line for the sake of simplicity. 


Generating Gar, r) 
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Economizing on the pre-processing stage 
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Matrix 
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Loved e The absorber can be hollow, 
BEMA | making the pre-processing 
Conditioning problem almost a surface one. 
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e A high degree of symmetry also helps. 
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e A given G,44 can be used for many problems. 
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Rule of thumb for avoiding the computation of annulled 
elements 







Two lines-of-sight (LOSs) are shown for a basis function at the 
point Py. These lines of sight skim the circumference of the 
(flattened, red) absorber. It turns out that only the interactions 
over the arc (P, : P : P5) between the two LOSs are needed as 
elements in the column representing the basis function at Pp. 
This results in substantial saving in matrix fill-up time. 
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Current distributions on conditioned and ill-conditioned 


circular cylinder problem 
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Symptom: condition number 
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JORE The graph below shows the condition numbers for the MoM 


Electric- 
Wende matrices: free space Green's function is red, BEMA is blue. 
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Green's Resonant frequencies, where the condition number blows up, are 
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@ Example: PEC Circular Cylinder with Diameter d = 3.1, 
Excited by a TM Plane Wave 
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Sanity check: using fully populated Gmod 
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Detail: The behavior of Gmoa (fully populated) for the case of 
the flattened absorber 





Detail: Five typical rows in the MoM matrix representing 
δ 


Balanced 
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Thinned using the line-of-sight rule - large circular filler 
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The number of non-zero elements in the matrix is then cut by 
about half along each row and each column, therefore the total 
number of remaining elements is about 2576 of the original 
matrix. Resultant current distribution is shown in the figure (blue 
line) compared with the conventional MoM solution (red line). 
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@ Conclusions 








A surface scattering problem can be solved using about 25% of 
the MoM matrix elements with little effect on the accuracy. 
Conditioning of the EFIE is also achieved. 


ə By inspecting lines of sight between basis functions it may 
be possible to avoid the computation of the matrix elements 
that represent severe attenuation. 

@ Solution is robust. 

e The effort incurred in the pre-processing stage can be 
alleviated as described above. 

e Most annulled elements need not be computed at all. 

The thinned matrix then forms the basis for either a direct 


method or pre-conditioned iterative method for 
formulating and solving surface integral equations. 
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